Signatures of resonance superfluidity in a quantum Fermi gas 
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In this letter, we predict a direct and observable signature of the superfluid phase in a quantum 
Fermi gas, in a temperature regime already accessible in current experiments. We apply the theory of 
resonance superfluidity to a gas confined in a harmonic potential and demonstrate that a significant 
increase in density will be observed in the vicinity of the trap center. 
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Following the successful realization of Bose-Einstein 
condensation (BEC) in confined vapors [Q, it is natu- 
ral to consider possibilities for observing the analogous 
superfluid phase transition in a dilute Fermi gas. Quan- 
tum degeneracy has already been demonstrated in a two- 
component Fermi gas of "^^K atoms ||^, although the spin 
states utilized were not suitable for exhibiting superflu- 
idity. The lowest temperatures achieved to date in this 
system are around 0.2 Tp — limited by Pauli blocking as 
well as a number of technical considerations j^] . In other 
experiments, the rethermalization of fermion atoms by 
elastic collisions with a bath of ultracold bosons is ex- 
ploited; realized by mixtures of ^Li and ^Li at Rice 
and at ENS j^, and more recently by mixtures of 
atoms and *^Rb atoms at JILA 

In order to observe a superfluid phase transition at 
critical temperatures as high as 0.2 Tp, the existence of a 
strong coupling mechanism which could lead to a sig- 
nificant amount of Cooper pairing is necessary. Sev- 
eral theoretical papers have presented models to inves- 
tigate this regime, essentially based on application of the 
Bardeen-Cooper-Schrieffer (BCS) [Q theory of supercon- 
ductivity. These approaches consider dilute Fermi vapors 
in which the two body scattering processes are charac- 
terized by a large negative scattering length a Q. Un- 
der such conditions the relevant length scale — the spa- 
tial extent of the Cooper pair — may become compara- 
ble to the average interparticle spacing. This places the 
system in a crossover region from the BCS superfluid- 
ity of momentum-correlated fermion pairs to the BEC of 
tightly bound composite bosons. In this crossover regime, 
fluctuations play a crucial role [|| and must be addressed. 

Eventually, as the coupling is increased, it becomes 
necessary to construct a theory in which explicit treat- 
ment of the composite bosonic states is made. Such 
an approach was proposed in the context of high- 
temperature superconductivity [lO| and is based on an 
effective many-body Hamiltonian, in which quasibound 
pairs are explicitly treated as resonance states embed- 
ded in the continuum of the Fermi sea. Such resonances 
are ubiquitous in atomic physics, where, for example, a 
Feshbach resonance |Q can be utilized to tune a qua- 
sibound state through threshold, providing an explicit 
microscopic basis for a theory of resonance superfluid- 



ity 

A convincing method for detecting the superfluidity 
will be required. Various approaches have been proposed; 
including measurements of the pair distribution |14f| , ex- 
periments involving the breakup of the Cooper pairs , 
measurements of the moment of inertia |16|] , and probes 
of collective excitations In this letter, we show 

that a more straightforward and direct experimental sig- 
nature of the transition to the superfluid phase is pro- 
vided by the density characteristics in an inhomogeneous 
system. We demonstrate that in a harmonic trap, the su- 
perfluid state is manifest as the appearance of a bulge in 
the central atomic density. To this aim we derive a the- 
ory of resonance superfluidity including the description 
of external conflnement. 
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FIG. 1. Real (solid line) and imaginary (dashed line) com- 
ponents of the T-matrix for collisions of the lowest two spin 
states of ^''K at a detuning of 20 Ef, shown in length dimen- 
sions, i.e. Tfe/(47rfi'^/m). The scattering length is the intercept 
at zero scattering energy which for this case is approximately 
-10000 ao, where ao is the Bohr radius. The large variation 
in the T-matrix over the relevant energy range indicates that 
a quantum field theory developed from this microscopic ba- 
sis will in general need to account for physics beyond the 
scattering length approximation. The inset shows the scat- 
tering length as a function of detuning, with 20 -Ef detun- 
ing indicated by the dashed-dot line. This curve obeys the 
following form a = abg(l — n/vo) where abg = 176 ao and 
K = 0.657 mK The quasipotentials to be renormalized 

are then Uo = 47rft^abg/m and go = \/nUo- 
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Although our results will have general applicability, for 
the purpose of illustration, we consider a typical system 
of = 5 X 10^ atoms in an isotropic harmonic trap with 
angular frequency u> = 27r 100 s~^. This gives a Fermi 
energy oiEp = {3Ny/^Huj. We consider a Feshbach 
resonance, illustrated in Fig. 1, for s-wave scattering of 
atoms in the lowest two hyperfine spin states which we 
denote by a G {Tii}- We begin by considering the gen- 
eral structure of the theory for the homogeneous system 
where the fermions are represented by the wave-number- 
k dependent annihilation operators a^^ and the compos- 
ite boson field by bk- The Hamiltonian is 



ha h 

+ ?7 ^ '^g/2+fcT^l/2-fci%/2-fc'i%/2+/c'T 
qkh' 



The self-energy has hermitian structure 



(1) 



kq 



where h.c. denotes the hermitian conjugate. The free 
dispersion relation for the fermions is e^, and v denotes 
the detuning of the boson resonance state from the zero 
edge of the collision continuum. The coUisional inter- 
actions include both background fermion scattering and 
interconversion between composite bosons and fermion 
pairs. It is implicit in treating U and g as constants that 
the theory will be renormalized, and thereby contain no 
ultraviolet divergences in the calculation of observable 
quantities po|] . This procedure involves ascribing a cut- 
off value K as the upper limit of all momentum summa- 
tions, and renormalizing the Hamiltonian constants in 
terms of K and the parameters for the Feshbach reso- 
nance as given in Fig. 1. Defining a — mK / {2Ti'^fi^) and 
a dimensionless factor F = (1 — aUo)^^ , the rcnormal- 
ization is executed by the following relations U — FJ/q, 
g = Fgo, and v ^ + cnggo- All results presented here 
have been shown to be independent of K . 

From this Hamiltonian, we construct the dynamical 
Hartree-Fock-Bogoliubov equations for both the bosonic 
and fermionic mean- fields. These equations involve the 
mean fields corresponding to the spin density n ~ 
'^k('^ka^ka) (ta-kcu to bc identical for both spins), the 
pairing field p = X]fc('^-fej.'^fc-f)i ^^'^ condensed boson 
field (j)m = {bk=o)- The single particle density matrix 



/ «fcT \ 



A = 



'hi 



(2) 



evolves according to the Bogoliubov self-energy S 

dg 



ih- 



dt 



(3) 



S = 



f Uk A 

Uk -A 

-A* -Uk 

V A* -Uk J 



(4) 



where the single particle energy is Uk — ek — n + Tkn, 
the gap is A = Up + g (pm , and ii is the chemical poten- 
tial. Here, without cost, we have upgraded the mean-field 
contribution to the single particle energy (which would 
otherwise be given by Un) to the full ladder sum, T^n, 
where Tk is the two-body T-matrix. This expression for 
the mean-field contribution to the particle energy is an 
approximation which is accurate for dilute Fermi gases 
where the quantum Fermi pressure limits achievable den- 
sities. The dynamical equations are closed by the evolu- 
tion equation for the boson mode; 



ih 



d(j),n 

dt 



9P- 



(5) 



The self-energy S is diagonalized locally at each k 
by the Bogoliubov transformation generating quasi- 
particles with energy spectrum Ek — \/U'k + |Ap. 
In equilibrium, the quasiparticle states are occu- 
pied according to the Fermi-Dirac distribution Uk = 
[exp((£'fe — ^)/fc6r) + l] . The corresponding max- 
imum entropy solution for the molecule amplitude is 
found by ihd(j)rn/dt = ^rn4'm where — 2//, so that 
Eq. (|^) implies = gp/{^rn — v)- The mean fields 
can then be determined by integration of the equilibrium 
single particle density matrix elements, given by; 



(27r 



1 

/ dk \{2nk - l)cos26lfe + l] 
■^Y Jo 



P 



(2^)^ Jo 



K 



dk {2nk - l)sin26lfe 



(6) 



where tan 29k = | A|/[/fe is the Bogoliubov transformation 
angle. Since 9k depends on n and p, these equations 
require self-consistent solutions 

So far this theory has been presented for a homoge- 
neous system, while we are interested in a gas of A^ atoms 
confined in an external trapping potential V{r). How- 
ever, a full quantum mechanical treatment of the trap- 
ping states is not required. For instance, in our case, 
with a temperature of T = 0.2 Tp, the harmonic oscil- 
lator level spacing is smaller than both the Fermi and 
thermal energies. Under these conditions, we may incor- 
porate the effect of the trap through a semiclassical local 
density approximation p2| . This involves replacing the 
chemical potential by a local one /i(r) — n — V{r), and 
determining the thermodynamic solution at each point 
in space as for the homogeneous system. 

In general, the validity of the semiclassical approxima- 
tion requires a slow variation in the occupation of the 
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discrete quantum levels as a function of energy. Remark- 
ably, in both bosonic and fermionic gases, this condition 
can often be satisfied even at very low temperatures; be- 
cause of strong correlations in a BEC due to repulsive in- 
teractions, and because of exchange effects in a quantum 
Fermi gas. In both cases the zero-temperature semiclas- 
sical approximation for dilute gases is usually referred to 
as the Thomas-Fermi approximation. 




Radial Position (nm) 

FIG. 2. Density profile at temperature T — 0.2 Tf and 
detuning i/q = 20 Ep showing accumulation of atoms at the 
trap center (solid line) . We compare with the profile resulting 
from the same fi but artificially setting the pairing field p to 
zero so that no superfiuid is present (dashed line). 

We evaluate the thermodynamic quantities at given T 
and N in three steps: (i) For given /i, we determine the 
local chemical potential /i(r) = ji — V{r) and use this 
value to find the self-consistent solution for the density 
n(r) and pairing field p{r) at each point in space, ac- 
cording to the solution of Eqns. (||); (ii) We modify the 
global chemical potential /i until the density integral is 
the desired atom number, i.e. N = J d^rn{r); (iii) We 
use the resulting solution for fi to calculate observable 
quantities, such as the density, gap, compressibility, and 
so forth. 

The resulting solution for the density distribution is il- 
lustrated in Fig. ||. A striking signature of the resonance 
superfluidity is evident in the predicted density profile 
which has a notable bulge in the region of the center of 
the trap. Experimentally, this signature appears to be di- 
rectly accessible. A typical approach would be to fit the 
expected density profile for a quantum degenerate gas 
with no superfiuid phase to the wings of the distribution 
(outside the dotted lines shown in Fig.^). Then the ex- 
cess density observed at the trap center can be recorded. 
Fig.^ illustrates the emergence of the superfiuid as the 
temperature is decreased. Qualitatively this situation is 
reminiscent of the central condensate peak observed for a 
Bose-Einstein condensed gas in a harmonic potential, al- 
though the connection appears to be somewhat serendip- 



itous. 




FIG. 3. Emergence of the coherent superfiuid for 
vo ~ "20 Ef- The superfiuid occupies an increasing volume as 
the temperature is reduced. Shown is the excess density (dif- 
ference between the dashed and solid lines in Fig.^) at each 
temperature. 

We explain the observed behavior by considering the 
compressibility of the normal and superfiuid gas. Ther- 
modynamically, the isothermal compressibility C is de- 
fined as C'-^ = n{dP/dn)T where P is the pressure, and 
is shown in Fig. ^ The compressibility is positive ev- 
erywhere, indicating that, in spite of the large attractive 
interactions, the Fermi pressure makes the configuration 
mechanically stable. A significant feature is the discon- 
tinuous behavior at the radius from the trap center at 
which the superfiuid changes from a zero to a non-zero 
value. This discontinuity is a manifestation of a second- 
order phase transition occurring in space. The disconti- 
nuity is a consequence of the local density approximation, 
and cannot occur in a finite system. However, a rapid 
change in the compressibility is expected. In principle 
this could be probed by studies of shock waves generated 
by the abrupt jump in the speed of sound as a density 
fiuctuation passes through this region. 

In conclusion, wc demonstrated that there exists a di- 
rect signature of superfiuidity in trapped fermion gases. 
The onset of superfiuidity leads to a density bulge in the 
center of the trap which can be detected by direct absorp- 
tion imaging. The critical conditions for superfiuidity are 
satisfied initially in the trap center, and the region of 
non-zero pairing field spreads out from the center as the 
temperature is lowered further. The increase in the den- 
sity profile in the superfiuid region is caused by a jump in 
the compressibility. Direct measures of this behavior are 
possible by the study of the propagation of sound waves. 
We have applied our method here to ''"K, but a similar 
approach is easy to derive for other interesting atoms, 
including in particular ^Li which is the other fermionic 
alkali currently being investigated experimentally. 
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FIG. 4. Inverse isothermal compressibility C^^ in units of 
the Fermi energy (solid line). Here fo — 20 Ef and T = 0.2 Tf 
(as can be seen from the limiting behavior at large radial 
position). A discontinuity appears at the radius at which the 
superfluid emerges (dotted line). We compare this solution 
to that corresponding to zero pairing field and no superfluid 
phase transition (dashed line). 
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